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Motivated by the well known impact of rare decays of hadrons and leptons on the evolution of the Standard 
Model and on limits for new physics, as well as by the proposal for Giga-Z option at TESLA, we investigate the 
rare decay Z —> bs in various extensions of the Standard Model. 


1. INTRODUCTION 

The central role played by rare decays on our 
understanding of elementary particle physics, is 
well known, where “rare” stands here for Flavor 
Changing Neutral Currents (FCNC), which are 
either small or practically vanishing in the SM. 
Some highlights: are: 

1. In K physics: The first appearance of charm 
in loops from which m c « 1.5 GeV was predicted 
[ 2 ]- 

2. In B physics: The importance of b —> sy in 
the SM and for extracting limits on Beyond the 
SM (BSM) scenarios [3]. 

3. The top quark FCNC provide an excellent tool 
to investigate various extensions of the SM [4]. 

4. The experimental upper limit of the decay fi —> 
ey [5], places severe limits on extensions of the 
SM. 

In view of the above and prompted by the re¬ 
cent discussion of a Giga -Z option at TESLA [6] 
in which the center-of-mass energy will be low¬ 
ered to M z , producing more than 10 9 Z bosons 
(i.e. ~ 100 times the number produced at LEP), 
one should investigate rare decays of Z s. Now 
since the important subject of rare leptonic Z —> 
hij, ii ij decays, for which the SM branch¬ 
ing ratio is < 10 -54 , was covered by Illana [7], 
we concentrate here on purely hadronic FCNC 
Z —> didj, di 7 ^ dj decays. In fact we only 
discuss Z —» bs which in most models, including 
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the SM, has the largest branching ratio among 
hadronic FCNC Z decays. Note however that 
experimentally, the latter is practically insepa¬ 
rable from the bd mode. Let us also note that 
when referring to the bs mode, we actually mean 
Z —> bs + bs. We note here that in the SM [8] 
BR (Z -> bs) « 3 x 1CT 8 . 

In the following sections we will discuss two 
variants of 2 Higgs Doublet Models (2HDM) and 
two of Supersymmetry (SUSY). Of the latter the 
first one will be: SUSY with squark mixing, while 
in the second one FCNC will result from SUSY 
with R Parity Violation (denoted by RPV, or 
j? p ). As we will see, BR (Z, —> bs) can be either 
smaller, the same or above the SM with a maxi¬ 
mal value of BR (Z —» bs) ~ 10 -6 . 

Experimentally, the attention devoted to 
FCNC in hadronic Z decays at LEP and SLD 
has been close to nil. The best upper limit is [9] 
£ 9 =d, s BR ( z -*■ b q) < 1-8 x 1(T 3 @ 90% CL. 
This is a preliminary DELPHI limit (which will 
probably remain as such forever...) based on 
about 3.5 x 10 6 hadronic decays. Experimental¬ 
ists who are privy to LEP and SLD data should 
be encouraged to look in their data and improve 
the above limit. 

Due to space limitations, the following discus¬ 
sion of various BSM models and their predictions 
for Br (Z —> bs), will be sketchy. Many more de¬ 
tails and a more complete set of references can be 
found in [1], In fact, almost each reference should 
start with: “See e.g ....” and end with: “... and 
references therein.” 




Figure 1. One-loop diagrams that contribute to 
the flavor changing transition V —> didj, due to 
scalar-fermion exchanges. 


2. GENERIC CALCULATION 

We start with a generic calculation of the di¬ 
agrams which modify (at one loop) the Vdidj 
vertex, due to charged or neutral scalar, as de¬ 
picted in Fig. 1. In our case V = Z, di = b 
and dj = s. The indices i.j and a, /3 indicate 
which fermions and scalars we are considering, 
respectively. The Feynman rules are: 

VJlfj '■ *7M ( < a L(Vf)L a R(Vf )R) 

V^SaSp : ig^f (p a - pp)^ 

S a djfi : i (l?i (a) L + 4 (a) R), 

where L(R) = [1 — (+) 7 5 /] 2 . 

There are 4 one-loop amplitudes, each corre¬ 
sponding to one of the 4 one-loop diagrams. 
Each amplitude is proportional to e^{q) times 


u(pb) [ 7 m {A l L + A r R) + ( B l L + B r R) p^} v(p s ). 
Al.Ri B r r are momentum dependent form fac¬ 
tors, calculable from the diagrams. There are 4 
per diagram, thus we have 16 form factors. Al 
for diagram ( 1 ) is: 

Al = 

and similarly for the other 15 form factors. C 24 
is one of the usual one-loop scalar functions [ 10 ] 
at m). , m 2 Sa , m 2 Sfj , m 2 di , q 2 , m 2 dj . 

Finally: 
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where A-[ is the Total sum of A R s from the 4 
diagrams, and similarly for A^, Bj^ and B 


3. MODELS AND PREDICTIONS 

The stage is now ready for identifying, for each 
model, the relevant scalars S a , fermions /j and 
the couplings a, b and g (with the appropriate in¬ 
dices), as expressed in the Feynman rules above. 
Then, the route for obtaining F(Z —> bs) using 
the generic equation in the previous section is 
clear. 


3.1. Two Higgs doublet models 

In 2HDM with flavor diagonal couplings of the 
neutral Higgs to down-quarks, the FCNC Z —> bs 
go through the one-loop diagrams in Fig. 1. The 
scalars are the charged Higgs bosons, S a —1 = H + 
and the fermions are /) = Ui, i = 1,2,3. The 
couplings are: Z^UiUj is as in the SM (therefore 
only i = j survives), Z^H + H~ is derived from 
the kinetic energy part of the Lagrangian C and 
H + iiidj is obtained from the Yukawa part which, 
in common notation is [ 11 ]: 


Cy = — 


^Qi[(up 1+ [/p 2 )4 


(Dl^i + D 2 $ 2 ) 4 


A choice of U and D, which are 3x3 matrices in 
flavor space, leads to a specific 2HDM. We now 
study two variants of 2HDM. 



3.1.1. Two Higgs doublet model of type II 

In this model, called 2HDMII, U 1 = D 1 = 0, 
Z —► bs was considered before [12]. Using real¬ 
istic values in the tan (3 — m R + plane, we ob¬ 
tain: BR(Z —■> bs) 10 -10 , two orders of magni¬ 
tude below the SM. 

3.1.2. Two Higgs doublet model “for top” 

In this variant [13], named T2HDM, the top 
is rewarded for its “fatness” by having its mass 
proportional to the large i> 2 , while all other 
masses are proportional to v \. It therefore makes 
sense to consider here only tan/3 >> 1. Using 
T2HDM parameters consistent with data we find: 
BR(Z -► bs) < 10 -8 . 

3.2. Supersymmetry with squark mixing 

FCNC in SUSY can emanate from squark mix¬ 
ing in: 


Cr" = -QKMq)hQj - U}(Mir)ijUj 

- Dl(Mp)ijDj + A»QiH u Uj + A i iQ i H d D j , 


with the usual notation for the squark fields [1] 
and where i,j are generation indices. Further¬ 
more, 
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certain assumptions [14] and taking only b — s or 
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The above Ss represent squark mixing from 
non-diagonal bilinears in C. mo is a common 
squarks mass scale, obeying: mo » Mz■ Also, 
Slrs will stand for squark mixing from non¬ 
diagonal trilinears in C [1]. For them we adopt the 
Ansatz of [15], leading to Slr oc vA/m%, where 
A is a common trilinear soft breaking parame¬ 
ter for both up and down squarks. Mfj v be¬ 
come 4x4 matrices in the weak bases $2) u = 
( s L ,s R ,b L ,b R ), (c L ,c R ,t L ,t R ). They are diago¬ 
nalized to obtain the mass eigenstates $d,u = 


(h,82,biM, (ci,C2,ti,t 2 ), with the help of 

Ru,d which rotates $ to <f>°. 

We can now describe two cases of squark mix¬ 
ing: b — s and t — c mixing. 


3.2.1. b — s mixing 

The scalars here are S a = $D,a, a = 1, 2,3,4, 
since b—s admixture states run in the loops. The 
gluon is the only fermion in the loops, thus /, = 
g. The a couplings are 0, since Zgg = 0. In 
other words, one diagram (out of the four generic 
diagrams) vanishes. The b and g couplings are 
functions of elements of the rotation matrix R r 
mentioned above [16]. Since the two Slr ~ 10 -2 
[17], we neglect them. For the other four Ss we 
assume a common value, i.e. o L £ = d L l = 
S RR 3 * = Srr 2 ' 1 = > an d va iT 0 < <5-° < 1. 

The parameters needed for masses, mixing and 
r(Z —> bs) are: mo, /u, A, tan/3, mg and 5 D . 
We vary them subject to m sqU arks > 150 GeV and 
have plots of practically everything as a function 
of everything [1]. 

We find BR(Z —> bs) ^ 10 -6 , where the highest 
value is attained for mg and one mj. ss the EW 
scale, while mj., j ^ i ss few TeV. Such split¬ 
ting requires “heavy” SUSY mass scale with soft 
breaking parameters, which is consistent with the 
non-observability of SUSY particles so far. 


3.2.2. t — c mixing 

In this scenario the scalars are S a = $u,a, ol = 
1,2,3,4, similarly to the previous case, except 
for D —> U. Obviously, t — c admixture states 
run in the loops. The loop fermions are the two 
charginos /, = yy, i = 1,2, and all four generic di¬ 
agrams contribute to Z —> bs. The Feynman rules 
[16] involve elements of the rotation matrix Ru 
mentioned above and the chargino mixing matri¬ 
ces. 

At the end of the day, running with the pa¬ 
rameters over all values consistent with the data, 
and with m,j > 150 GeV and m x > 100 GeV 
we obtain: BR(Z —» bs) 4$ 10 -8 , which we could 
have anticipated since BR(through t — c mixing): 
BR(through b — s mixing) ss a : a s . 



3.3. Supersymmetry with RPV 

Since there is no sacred principle which guar¬ 
antees R-parity conservation, we assume in this 
part of the talk that Rp is violated. Then, l/l p 
terms in the SUSY superpotential W lead to 
FCNC. A terms (pure ff) in W are irrelevant 
at the 1-loop level. In addition we assume, for 
the pure /B terms, that A" << A', and also 
neglect the bilinear term in the Tft p part of W. 
Then: Wrpv = ^abKjk^iQ’j^k- In addition, if 
1/tp is OK then the Rp conserving soft SUSY 
breaking is extended. We need only the bilinear: 
Vrpv = eabbiL-Hl where L, H u are the scalar 
components of the hatted L and H u , respectively. 
We therefore have two types of FCNC: 

Type A : Trilinear-trilinear: F [Z —> bs) oc 

(A'A') 2 . 

Type B: Trilinear-bilinear: T(Z —> bs) oc (b\') 2 . 

3.3.1. Type A: Trilinear-trilinear terms 

We further sub-divide type A contributions to 
6 groups according to the scalars and fermions 
running in the loops. For instance, in type A1 
the scalars are S a = eL )0n a = 1,2,3 and the 
fermions are /,; = iq, i = 1,2,3. The a cou¬ 
plings are identical to their SM values, 6 l = 0 
(for all i,j and a), b % g [a) = — A£Y and g’f = 
—e(cfr — s^ v )8 ol p/2swCw- Unfortunately, going 
over all type A groups, taking into account the 
available limits on A's and on the other relevant 
parameters, we obtain for the trilinear-trilinear 
case: BR(Z —> bs) 10 -10 . Our results are in 
agreement with the special cases in [18]. 

3.3.2. Type B: Trilinear-bilinear terms 

In this case, a Higgs exchanged in the loop 
mixes with a slepton, through e a bb 3 L 3 Hu, as¬ 
suming that only b 3 ^ 0. We choose to 
work in the “no VEV” basis v 3 = 0 in 
which: H u = (h+, + v u + i<f>°)/V 2), H d = 

(KS + ^ + ^i/V 2 , h d ), 

L 3 = ((£ 9 + iv°_ )/v / 2 , eg ) , where v p , i>°, e 3 
are SU(2) CP-even, CP-ocld r-sneutrinos, fp, re¬ 
spectively. In the basis = (h+ ,h^,e 3 ) we 
wrote the mass matrix in the charged scalar sec¬ 
tor, in the basis $° E = (£°, P+) we wrote the 

mass matrix in the CP-even neutral scalar sector, 
and in the basis <J>q = (</> 9 , </>°, ) we wrote the 


mass matrix in the CP-odd neutral scalar sector. 

The new charged scalar and CP-even and CP- 
odd neutral scalar mass-eigenstates are obtained 
by diagonalizing the above-mentioned matrices. 
They are: = {H + , G + , f+), = ( H,h,v £), 

and <1>o = (A, G°,vZ.)- In the limit b 3 —» 
0: H , h, A , H + become the usual ones. Rotat¬ 
ing with the diagonalizing Rc,e,o (for Charged, 
Even-CP, Odd-CP) matrices, one goes from the 
$s to the <J>°s. All depends on the four parame¬ 
ters A 0 , m° sv (the masses in the limit b 3 —> 0,), 
b 3 and tan/3. 

Let us sub-divide type B into two types accord¬ 
ing to the scalar and fermion in the loop: 

Type Bl: Here S a = 'Lc.cd fi = u i with a = 
1,3; i = 1,2,3. The a couplings are equal to 
their values in the SM. The b couplings include 
elements of the rotation matrix (for the charged 
fields) Rc and A', and g 1 ^ 3 = —ecot2 Owba/3- 
Type B2: In this case S a = and <j>o,/ 3\ fi = 
di with a = l,2, 3, and /3 = 1,3; i = 1,2,3. This 
is the only case for which our generic form is insuf¬ 
ficient. This fact results from the appearance of 
two new diagrams proportional to a scalar-vector- 
vector coupling (ZZ&e in our case). The other 
eight diagrams are special cases of the generic 
ones in Fig. 1. 

Inserting parameters consistent with the data 
we found that for type B: Br (Z —> bs) F, 10 -6 . 

4. EXPERIMENTAL FEASIBILITY 

Let us briefly comment about the feasibility 
of observing (or limiting) a signal of BR(Z —> 
bs) ~ 10” 6 , at a Linear Collider producing 10 9 
Z-bosons. Such a signal leads to one 6-jet and 
one q- jet, where q stands for quarks lighter than 
b. The main background is from Z —> bb. Us¬ 
ing what, we think, are realistic efficiencies we 
find that a new physics signal Z —> bs, with a 
branching ratio of order 1CU 6 , can reach beyond 
the 3-sigma level [1]. We can also get a clue about 
how low one can go in the value (or limit) of 
BR(Z —* bs) with 10 9 Z-bosons, from the fact 
that the DELPHI preliminary result reached [9] 
BR{Z —> bs) < e>(10 -3 ) with 0(1O 6 ) Z-bosons. 
Scaling this limit, especially with the expected 
advances in 6-tagging and identification of non-6 



jets methods, an 0(1O~ 6 ) branching ratios should 
be easily attained at a Giga -Z factory. One needs 
realistic simulations as feasibility studies for this 
important rare Z decay mode. 

5. SUMMARY AND CONCLUSIONS 

Our results are best summarized in Table 1 
which shows the best values for Br (Z —> bs) in ex¬ 
tensions of the SM discussed in this talk. The SM 
result is given for comparison. Note that we have 
not included interference with the SM as each of 
the values “stands alone”. In some cases such 
interference may increase the branching ratio to 
~ lCT 7 . 


Table 1 

Maximal By(Z —> bs) in extensions of the SM 


model 

scalars in loop 

Br 

SM 

W + (no scalars) 

3 x 10 -8 

2HDMII 

H+ 

10- 10 

T2HDM 

H+ 

10" 8 

b-s mix 

b-s admix 

10 -6 

t-c mix 

t-c admix 

10" 8 

tri-tri 

q, d,£ 

10- 10 

tri-bi 

v-h, H 1 A admix 

10 -6 


We conclude that Giga -Z experiments will have 
the opportunity to place significant limits, or 
hopefully discover the scenario beyond the SM, 
by searching for hadronic (and leptonic [7]) neu¬ 
tral current flavor changing transitions. 
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